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1. Introduction. 



The symmetry between low and high temperatures in finite temperature field 
theory has re-emerged recently in connection with the questions of entropy bounds 
and the Verlinde-Cardy relation. Some comments were made in an earlier paper 
[1] which we wish to enlarge on here. Although the main point of this earlier 
communication was the thermodynamic significance of any zero modes, this will 
occupy us only incidentally here. 

The primary emphasis in the present article is on inversion symmetries in finite 
temperature field theory and its relation to elliptic functions. Inversion symmetry 
usually appears as a consequence of Jacobi's inversion identity for theta functions, 
in one of its many guises, and it should be no surprise that elliptic functions feature 
in our discussion. 

In Section 2 we discuss for which manifolds simple inversion properties might be 
expected and we are naturally led to the manifold S^ x S'', S^ being the thermal circle 
and S'^ the d-dimcnsional sphere. A unified approach to the energy of quantum fields 
in curved space-time is via the thermal ^-function and some basics are presented 
in Section 3. Afterwards we relate the energy to Eisenstein series which allow 
symmetry properties to be discussed very easily. In particular, we can use relations 
between the Eisenstein series and the Weierstrass p-function to find the scalar 
energy for all dimensions just in terms of the results on the three and five sphere. For 
spinors one iterates from the one and three dimensional spheres. Elliptic functions 
are then used in Section 6 and 7 to give a description of temperature inversion in 
terms of the modulus, k. The relation to elliptic functions also allows one to express 
the energy in terms of powers of k'^ and of the complete elliptic integral, K{k). This 
has the advantage that, at temperatures corresponding to singular moduli, closed 
forms for the energies can be found, see Section 9. The Conclusion describes the 
main points raised in the article. Three appendices contain technicalities about 
different viewpoints on the approaches for spinor fields, power series expansions of 
elliptic functions, and some explicit expressions too cumbersome for the main body 
of our paper. 

2. Background technical facts. 

Any symmetry relation has to be expressed in a dimensionless quantity say 
aT = a/ 13 = ^/27r, where a is some length scale associated with the spatial manifold, 
M.. Clearly reciprocal symmetry involves exchanging (5 and a. 
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Such a symmetry is therefore suggested whenever one can detect the torus 
structure, S^xS"^, in the thermal manifold, x M, the above exchange corre- 
sponding to a transposition of the two circles. An example where the symmetries 
are clearly visible is when Ai is the torus, S^xS^ x . . . xS^. However the case that 
interests us here is when A4 is the d-sphere, in particular d odd. The sphere, and 
all spaces of constant curvature, have only one relevant length scale, which at least 
simplifies the possibilities for reciprocation. 

A basic analytical fact is that the scalar spectral quantities such as the heat- 
kernel (or quantum propagator in our case) on the rf-sphere can be 'integrated' in 
steps of two down to that on the circle, for theories conformal on the space-time. 
The thermal space, S^xS'^, reduces to S^xS^. Alternatively, the thermal circle, 
S^, can be differentiated to give a "thermal d-sphere" resulting in the symmetrical 
combination S'^xS''. This occurs only if d is odd and explains why inversion sym- 
metries arise in this case. For even d, the scalar heat-kernel can be fractionally 
differentiated to give S'', but only if it is twisted, i.e. a ^2~function, e.g. [2]. This 
is clearly the origin of Cardy's use of anti-periodic conditions for bosons, [3], for, in 
order to achieve any symmetry, the thermal circle must be twisted as well. 

A practical strategy is to treat the terms in the degeneracy polynomial inde- 
pendently and deal with them one by one. Each part will have a different inversion 
property. It is this approach that will occupy us mostly in this paper. 

We might mention that these calculations can be extended to split-rank sym- 
metric spaces of which odd spheres are the simplest examples. (See, in this connec- 
tion, [2].) The heat-kernel expansion terminates at the firstterm on all (semi)-simple 
group manifolds, but only if the scalar is conformal in 4 dimensions, whatever the 
dimension of the group. 

3. Statistical mechanics. 

The essential points were made in [1] but have to be repeated here for com- 
pleteness and ease of exposition. A unified approach is via the thermal ^-function, 
C(s, /?), [4-6], which neatly incorporates any ultra-violet and zero mode effects. 

The general form of Ci^iP) is, on an ultrastatic space-time, TxM., 



where u;^ and d^ are the eigenvalues and degeneracies of minus the Laplacian, — A2, 
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(1) 



n 
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on the spatial section, M.. The dash means that the denominator should never be 
zero. 

The free energy (effective Lagrangian) is obtained from the limit 



F = - lim 
2 s^o 



(2) 



s 



which, in general, can have infinities. Leaving these aside, in [1] the question of 
inversion symmetry was discussed directly, and most conveniently, from (1) as dis- 
playing most clearly any symmetry under interchange of (3 and a scale a, which is 
embedded in the uj^, typically uj^ ~ 1/a. 

It is possible to write the formalism in several equivalent forms often charac- 
terised by the way in which the mode information is organised. Each form has its 
own features. In [1] two such were outlined and it is the one which starts from the 
familiar statistical sum that we wish to concentrate on here. 

In order to avoid certain problems, we deal only with internal energies, and 
leave for future work a treatment of the free energy and a complete thermodynamical 
scheme. 

As shown in [4] , the limit (2) leads, or can lead, to the expression for the total 
internal energy, E = d{PF)/dp, 



Eq is the zero temperature value and may be the result of a regularisation/ renor- 
malisation procedure. The effect of zero modes has also been included as the 
second term [7,8]. It is this form we wish to dwell on for a while as it allows one to 
make contact with some standard results in elliptic function theory. 

4. Scalars on the odd sphere. 

Information on modes on a sphere can be found in many places going back 
many years. For scalars conformal on the space-time, TxS*^, the eigenvalues are 
perfect squares and the information can be displayed in the (zero temperature) 
^-function. 




(3) 




(4) 
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where, as mentioned in the Introduction, we are focusing on odd spheres. The 
degeneracy is a Gegenbauer polynomial. Another way of writing this involves the 
Barnes ^-function but we will not need this here. 

Except for the circle, d — 1, there are no zero modes so for simplicity let's 
restrict initially to d > 1 and make the common expansion of the degeneracy (e.g. 
[9]), 

in' - (r - 1)^) ...n' = n'^-' (5) 

which allows ((s) to be written as a sum of Riemann ^-functions, as has been done 
many times before. Our object is not to investigate this particular aspect much 
further but only to say that from this sum, or otherwise, it is possible to obtain 
the finite values of the zero temperature vacuum energy, Eq, as a, sum of ordinary 
Bernoulli numbers, a calculation done a long time ago, and a number of times since. 
This is not the most compact form, but it is suitable for our purposes where, as has 
been indicated, we are going to treat the terms in (5) separately. For the {2t — 2)th 
power we have 

«^oW = ^Ci^(l-2t) = -^, 

according to general theory [4,10]. (There is no need to use the Bernoulli number. 
One could just leave the Riemann (^-function, which is actually better.) 

Let us therefore return to (3) and consider the 'partial' internal energies, 

aE(t) = -^^j:-^^, (6) 

n=l 

and, in order to produce an expression that looks familiar, define a new parameter, 
Q, by, 

q = e--/^ , (7) 
to rewrite the internal energy (6) as a q-senes, 

R ~ „2t-l„2n 



n=l 



The summation is an example of a Lambert series and occurs in analytic num- 
ber and in elliptic function theory. In the former, q usually is the square of the q 
here, (7). 
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The most rapid way of proceeding is to note the standard result that the 
combination (8) is, up to a factor, an Eisenstein series. Precisely, 

^t(o = (-irf^G,(i,z/o, (9) 

where the Eisenstein series is defined by (we are following Hurwitz, [11,12], here), 

oo 

Gti^^.uj^) = 7 ^ ■ (10) 

= — oo 

When t = 1, the torus case, the sums have to be considered more carefully. 
This point will be returned to. A useful discussion is given by Rademacher [13]. 

The behaviour of e^(^) under inversion is now apparent from (9) and the struc- 
ture of the Eisenstein series (10), i.e. 

et(i/o = (-ir^et(e). (11) 

The Gj. are holomorphic modular forms, invariant under SL(2,Z) action on the 
periods a;i,a;2. The simple inversion case ^ 1/^ produces (11). The more general 
action contains other information. 

The result (11) when re-expressed in terms of the expression in (8) is equivalent 
to the duality relation, 

"* E il;r3T -(-"')' Eil^ = ("'-(-"')') ^"'' = ''^ (12) 

n=l n=l 

which is given by Ramanujan, typically without proof. It was derived later by Rao 
and Ayyar, [14] and thereafter by many workers including Malurkar [15] and Hardy 
[16] who use Mellin transforms. Berndt [17] gives some history which we could 
modify by saying that, since (12) follows immediately from (9), then it is contained 
in Glaisher, [18] and, by the same token, appears even earlier in Hurwitz [11,12]. 

The Eisenstein series enter into the Laurent expansion of the Weierstrass p- 
function with periods uji and UJ2, 

^ 00 

p{u)^-+J2('^t-l)Gy-\ (13) 

t=2 

which can be regarded as a generating function for the internal energies, e^. 
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Because of the restrictions of analyticity and periodicity, p satisfies the Weier- 
strass differential equation, (originally obtained by Eisenstein), 

= 4p{uf - g2{u;^,uj2)p{u) - g^{u;^,uj2) 

which leads to a recursion formula giving all the as algebraic polynomials in just 
the two invariants, g2 and g^, with rational coefficients, which is a very remarkable 
fact. A few examples are given in [12], footnote on p. 33 and a larger list in [19] IV 
p.88. 

The calculation has been kindly done for us by Rademacher, [13], and gives, 
expressed in the e^. 

This is an identity in ^ and yields several specific identities generally of arithmetic 
interest. For example, in the low temperature limit, ^ = 0, i.e. q = 0, (14) reduces 
to an identity between Bernoulli numbers, e.g. [13] p. 124. 

The recursion is more simply expressed in terms of the expansion coefficients, 
= {2t - 1) G^, in (13), e.g. [20] p.92, [19] Vol I, p.l76, [21] p.l55, [22] p.ll, as 

3 



T.= -, -J2t,T,_i. (15) 

1=2 



(t-3)(2t + 



Examples of the minimal polynomials resulting from (14) are 

.... 2 2^3257 2^32527, 3 

€4=12062, €5 = — — — €263, €6 = — (96e2 + e3), 

08 34 r3 j2 

^7 = 2^335276163, 68 = 62 (44 6i+62). 



(16) 



The total energy on the d-sphere is, according to equations (3), (5) and (6), a 
sum of the e^, 

V 

a£;s.^:E,(0 = ^$^(-l)*+>t(Oet+i(0, d = 2r + l, (17) 

and therefore can also be expressed as a polynomial in 62 and 63, or, equivalently 
in the three- and five-sphere energies, and E^,. 
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This has numerical import since it means that the energies on all d-spheres are 
known once those on and are and this is true for all temperatures, including 
zero. For example, the energies up to the nine-sphere are, 

E3 = €2 
— 1 / 

_ 1 (18) 

Eg = (5040 £2 £3 - 18480 + 539 €3 - 396 e^) . 

Apart from this, the behaviour of the energy on under the inversion ^ — > 
1/^ is determined from (17) and the individual behaviours, (11). In this way we 
regain the high temperature form derived generally in [4] in terms of the heat-kernel 
coefficients. It is another way of introducing, organising and using the spectral 
information. 

Quantities equivalent to the are considered by Ramanujan [23] who computes 
the expressions analogous to (16) by a recursion derived purely algebraically. 

5. Spinors on the odd sphere. 

The spinor spectral data on the sphere are well known so we can be brief and 
simply display them via the zero-temperature (^-function on odd spheres, 

^ _ 2(.<+i)/2(_i)(.-i)/2 C;^/^^ ^,(n+ 1/2)^' 

where the degeneracy has again been expanded and some re-arrangement of the 
starting points made for convenience. 

We will again treat the individual powers in the degeneracy separately, and, 
putting some constants aside for ease, at finite temperature the fermion partial 
internal energies are defined by, 

^.(0 = -2C.(-l/2) + E eUnti)/^ + 1 ' (20) 

n=0 

where the first term is the zero temperature part with, 

Ct(s)=22*-2«-2^^(2^-2t + 2,l/2), 
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which can be given in terms of BernouUi numbers, if desired, as 

-^22*-'Ch(1 - 2t, 1/2) = (1 - 2^*-^) ^ = 22*-iMlM . (21) 

The point in writing it in this particular way is that we can now immediately 
appeal to a formula given by Glaisher, [18] p. 64, which reads, after transcribing the 
notation and reverting to the double-sided summations, 

{Zn+1) q ^ _(-\n2t-l\^ , 
l+g2n+l ^ ^ At 

n=0 

(-l)*+i(2t-l)! ^/ (-1)^1+^2 



'^'^'^^ mj^.m^ {rrii + 17712/$,) 



2t ' 



where q = e A glance at (20), (21) and (22) shows that the spinor partial 

energy is given by a sign-modulated, doubly twisted Eisenstein series, i?^, 

miO = -^^t — -Ht{i,i/0 (23) 

where, 

^t(^i,^2)= E (24) 

= — 00 

and obviously, therefore, enjoys the same inversion properties as the scalar quantity. 



^t(l/0 = (-l)*^^t(0- (25) 
Re-expressed by (20), this result becomes. 



n=0 n=0 

(26) 

for = TT^. According to Berndt, [24], this formula is due to Rao and Ayyar, [14], 
but we see that it is contained in Glaisher, [18] (but not in Hurwitz [11]). 

Without going into details, we note at this point that relations more general 
than (12) and (26) hold when t is allowed to be negative. 

As an already discussed example, [25] , take the three-sphere. The total internal 
spinor energy is, 

4(e) = k^2(e)-r7i(e)). 
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Under inversion, 



In this case the heat-kernel expansion terminates at the second term. 

The twisted Eisenstein series, H^, (24), appear as coefficients in the Laurent 
series for the twisted Weierstrass function, p, (an 'elliptic function of the second 
kind') defined by. 



oo 



This will appear later in Appendices A and B. 

However, in order to derive fermion equations analogous to the recursive (18), 
it seems that Jacobi functions are more suitable and so we turn to this description. 

6. Elliptic functions and alternative descriptions. 

In this, and the following section we expand on the previous analysis and also 
discuss inversion relations using different ingredients. 

Hurwitz proves (9) directly, starting from the partial fraction expansion for the 
cotangent. The proof in Rademacher [13] involves a slight reorganisation and an 
application of Lipschitz' formula, 

^ ' m=0 n= — oo ^ ' 

to the 777-1 summation. Equation (28) is a duality relation and can be proved using 
Poisson summation. There are some restrictions. Thus Re ,2>0, 0<q;<1 and 
Re s > 1. If < a < 1 then Re s > 0. 

The result, (9), is also given in Glaisher, [18] §72, and Halphen, [20], proved in 
slightly different manners. Halphen gives some specific examples on p. 446. (There 
is a minor error of 4 in the third line of his eqn.(78).) See also Hurwitz [12] pp. 22, 33. 

Halphen employs Weierstrass' whereas Glaisher develops everything from 
the Q'-series for Jacobi elliptic functions. The ones relevant here for bosons and 
fermions are those for zs and ds, the former being Glaisher's notation for the Jacobi 
(■-function, Z^. The following trigonometric expansions are standard and are due, 
effectively, to Jacobi; 

— zs « = cot — + 4 ^ ^ ' (29) 

n=l 
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and 

'IT/' 1 2n+l 

±!ldsu^ , , , - 4 y ^ , sin ((n + 1/2)^) , (30) 

71 Sin(7r«/2i^) ^l + g2n+l W ' ' R'' ^ ' 



n=0 

valid for \lmu\ < ttKqK' /2. 

The notation is the usual elliptic one, with, 



and K' given by, 
where /c^ + /c'^ = 1. 

In this scheme the key fact is that expansions of (29) and (30) in powers of 
u will produce expressions for the quantities and r]^, of (8) and (20), in terms 
of k and K. This is done by Glaisher who then proceeds to derive the relation to 
Eisenstein series, (9). Explicit forms are given later. For now we note the Laurent 
expansion of ds. 



^ (2/ - 1)! \2K 



exhibiting the partial energies, ry^, as coefficients. 

Temperature inversion is just the interchange, k k' so that K K' and the 
effect can be deduced directly from (29) and (30) because the functions zs and ds 
reproduce themselves, i.e. 

zs{u, k^) — —izs{—iUj k''^) — 



2KK' (33) 
ds(tt, k'^) = -i ds{-iu, k''^) . 

We might term this the Jacobian approach to inversion, and the one via Eisenstein 
series, the Weierstrassian approach. 

Some inversion formulae are given in Whittaker and Watson [26] p. 535 ex.60 
derived from the imaginary transformation of sd. See also Hancock [27], p. 308, 
Ex.4. 

The expansions (29) and (30) yield, on setting A; — > A;', w — > — m, 

zs(-z«, k'-") = coth ^ - 4 smh (34) 

n=l 
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and 



ds( -iu, k''^) = , , , 

it: ^ ' smh{'Ku/2K') 



r + 4y sinh((n+l/2)^), (35) 



so that from (33) and (34) we obtain the general inversion formula, 
K 



/2r7 

COth ^ - 4 2^ ^ _ ,2n 
n=l 



2K' 



K' 



cot 



TTtl 



2n 



72n 



sm 



n=l 



(36) 



+ u . 



which reproduces (12) on expansion in ti, after noting that ^ = nK'/K and //' = 
TiK/K' . We have said before that this is equivalent to the inversion symmetry (11). 

As a special case, an elementary consequence of (11) is that 6^(1) vanishes if t 
is odd, c/Glaisher [18], §§86-94. (This follows more directly from a consideration 
of the Eisenstein series.) It gives the interesting identity for the Bernoulli numbers. 



n=l 



B 



2t 



,2wn ]^ 



4t 



t = 3,5,... 



(37) 



rediscovered many times. Watson's proof 
sion for ds^ similar to (29). 

When t — 1 the identity is modified to 



involves an elliptic function expan- 



E 

n=l 



n 



1 _ 1 1 

Stt ~ 24 ~ Stt 



(38) 



The 'extra' contribution of l/87r arises from the inhomogeneous part of the 
imaginary transformation of the (^-function, zs, (33) and corresponds, in our earlier 
terminology, to the zero mode on the circle, S^. For this special case the vacuum 
energy is. 



B 



n 



n=l 



(39) 



The identity (38) can be derived via the full duality relation obtained from 



(36), 



B, 



/2n 



B, 



^Ei 

n=l 



nq 



2n 



,2n / + 27r ' 



(40) 
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by setting fi = fj,'. Equations (36) and (40) seem to have been derived first by 
Schlomilch, [29]. Curiously he does not produce (12) although it is implicit in his 
formula. 

Expressed in terms of the internal energy, (40) becomes, 

E,iO-^ = -eE,{l/0, (41) 

and we see that the zero modes have disappeared from the left hand side. 

In terms of the Eisenstein series, (9), one encounters the series Gi which for- 
mally diverges. Giving this term meaning, as described in Rademacher, for example, 
would produce the same result as above, but we have managed to bypass this analy- 
sis. Indeed, from (41), by recognising that the quantity possessing simple inversion 
properties is E-^^{^) — ^/47r, one can derive the expression for Gi, [13] (63.82), 

G.(M/« = (2.)^(f-2El^ + |)- 

^ n—1 ' 



7. Spinors again. 



In this section we discuss the corresponding equations for fcrmions and also 
give results analogous to the recursive formulae (16) and (18) for bosons. To begin 
with, we present the fermion equivalents of the inversions (36) and (37). 

The general inversion formula following from (30), (35) and (33), 



K 



sinh {T^ujlK') 



-,/2n+l 



n=0 



_|_ q'2n+ 



^ sinh ((n + 1/2) p) 



2n+l 



sin {ttu/2K^ 



-4Ey^^«-((^+V2)^) 



n=0 



is equivalent to the Fermi-Dirac relation, (25) or (26). The identity analogous to 
(37) is, 

(2n + l)4*+i ^ 1 ..S. 

g(2n+l)7r _|_ 2 4 ^ 



E 

n=0 



1) 



'4t+2 

2t+l 



(42) 



There is no zero mode special case. The absence is typical for spinor thermody- 
namical quantities in general. 

In order to extract the r]^ from (32), the explicit Laurent expansion of ds is 
required. A number of ways of obtaining this are described in Appendix B. The one 
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most appropriate at the moment depends on the differential equation satisfied by 
ds. The point is that the cubic recursion thence derived in Appendix B allows any 
r]f to be expressed as polynomials in iji and 772- This means that, in this case, the 
energies on the circle and thrcc-sphcre are sufficient to determine all energies. The 
reason for the contrast with scalars is the absence of an obstructing zero mode. 
The equations corresponding to (16) are, 

240 . „ 2\ 
40 

V4 = Y ('^^2 + ^^0V!V2 + 576174) (43) 
V5 = ^Vi (85 4 + 960 vh2 + 4032 Tjf) , 
and the spinor energies are, for example, 

— f 1 

^3(0 = 2(^2-^1), (44) 

4(0 = ^(9^-10% + %), 

f f f 

allowing to be expressed in terms of -E'i(0 -^slO, this may be continued. 



8. Explicit formulae. 

The value ^ = 1 is the lemniscate case and for scalars, e.g. [20] p. 64, 

6^(1) = r^(V4) e(i) = o. 

In elliptic function terminology, the lemniscate value corresponds to K = K' and 
so to a modulus of 1/ y/2. 

This is a specially interesting point, but it is an important fact in the practical 
application of elliptic functions that, given the value of K/K' , one can calculate the 
values of K and K' separately, and also of k. The physical significance of this for 
us is the following. 

Expansion of the right-hand sides of (29) and (30) in powers of u introduces 
the quantities and 77^, as has already been used to derive inversion symmetry. 
However now we see that the left-hand sides will yield expressions for and rj^ in 
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terms of k and K. These have been given by Glaisher, for example. Hence, knowing 
the temperature, i.e. ^ = K/K', one can compute the partial energies from elliptic 
function numerics. Furthermore, according to our previous analysis, see (18) and 
(44), only €2, €3, rji and 773 need to be evaluated in this way in order that all spheres 
are encompassed, in principle. 

Making use of Glaisher's expressions, for convenience, we find, for bosons, 

e, = -^-,{1 - + k') 

(45) 

= Q^C^ - k')i2e - + k') K\ 

As might be expected, the expansions occur elsewhere. The relation between 
the Eisenstein series and the elliptic functions follows from the very basic connection 
between Weierstrass and Jacobi, e.g. [20] p. 24, 

4K'^ 4 ^ ^ 

and power series expansion of the right-hand side again yields the {k, K) polyno- 
mials in (45). The expansion of is given by (13) as a power series in §2 and grg, 
the standard expressions for which are, [20], pp. 59-60, (with a sign error), [19] IV 
p.92, [30] p.152, [27] pp.201, 298, [31], p.57, [32] p.384, 

g,^^X\l-e + k') = ^X'{l-{kkr) 

93 = |a3(1 + P)(l - 2k^){l - k-/2) = Aa3(/,'2 _ J,2^^, ^ 



where A is defined by 



A= ( — 

UJ 



2K^ ^ 



Note that our periods are uj^ and 0^2 we are setting uji = 1 and U2 = iK'/K. 
For useful relational facts see [19] vol.11, chap.IV, pp. 180,207, and [20],pp.44-46. 

It is interesting to confirm the inversion symmetry directly from the forms in 
(45). Write generally e^{k,K) = F^{k'^) K"^^ (k) . Inversion is k"^ ^ k'^ = 1 - k"^. 
Then (Hancock [27] p.308, Ex.4), 

F,{1 - k') = i-iy F,{k') (48) 

so that, 

2t 



e,ik',K') = i-iy(^^^ e,{k,K), 
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which is exactly (11). 

Turning now to the fermion case, the expressions we need are most easily read 
off from the expansion for ds given in Glaisher, [18], 

% = (1 - 2e){16k^ - 16k^ + 31) , 

showing, as a simple check, that ?7i(l) = 0, ^(l) = at the lemniscate point, 
k = 1/ V2. The expressions satisfy the relation in (43) . 

The power series expansions of the elliptic functions used here, although old 
and standard, are not emphasised in the usual textbooks, with a few exceptions. 
For this reason they are developed a little in Appendix B, where a recursion for 
the expansion coefficients of ds is given and has been checked against Glaisher's 
expressions. 



9. Numerical evaluation. Singular moduli. 

Computationally, the nome, q, is known from the temperature so the modu- 
lus, k, can be calculated, usually via a power series in q, and then K[k) follows 
straightforwardly by the Legendre-Landen-Gauss method {e.g. Appell and Lacour, 
[32] §198, Hurwitz and Courant, [21] p.242, Greenhill [31], p.322) or one could, 
pragmatically, just use the built-in algorithms of Mathematica. 

Generally, there seems no advantage in using an elliptic formulation, as opposed 
to a simple evaluation of the defining thermal state summations, (8) and (20), taken 
in conjunction with inversion, (11) and (25). However, as well as the lemniscate 
case, there are many other special values of K'/K that have 'exact' expressions for 
k so that only the computation of K from k remains and in many cases, like the 
lemniscate one, K is known in closed form. 

Abel has shown, for K' / K of the special form, 

— = ———, a,h,c,d,neZ, (50) 
K c + d\/n 

that k is the solution of an algebraic equation with integer coefficients. Relatedly, 
it is known that this also applies when K'/K — \/iV, where is rational. In 
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particular, when N is an integer, the modulus, k = kp^, is referred to as a singular 
modulus. 

There are many classic works on modular equations and singular moduli. A 
little known one that gives a direct algebraic construction for prime order is Russell, 
[33]. Lists of fcjY can be found, e.g. in Weber [30,34] and in Greenhill [35]. The 
notebooks of Ramanujan also contain many examples of modular equations, [36]. 

It is also known that the complete elliptic integral, K, evaluated at a singular 
modulus can be expressed in terms of Gamma functions with rational arguments. 
The = 1, 3, 4 cases are classic, {e.g. [26]) and the N = 2 case is treated by Glasser 
and Wood [37], although it seems to be implicit in work of Ramanujan. 

The general case was discussed by Selberg and Chowla, [38], who proved that 
the complete elliptic integral, K{kj^f), is always expressible in terms of Gamma 
functions and gave the N = 5 and N = 7 forms. 

The lowest value. A?" = 1, is just the lemniscate value touched upon earlier. 
When N is a square, K{kj^) is an algebraic number times K{ki). More generally 
when N contains a square, K{kp^) is an algebraic number times K{k^^) where A^q 
is the square-free part of N. 

Much effort has gone into finding K{kj^) together with kj^ and, by now, many 
specific examples have been listed e.g. [37-45]. Borwein and Zucker, [44], show that 
using the beta function often provides simplified expressions. 

In the following we calculate closed forms of the energy in one, three and five 
dimensions for a selection of values, ^jv: of ^ that correspond to singular moduli. 

Examples we will look at are = 1, 2, 3, 6, 7, 10, 15, all of which can be found 
in Zucker and Joyce, [43] , for example. For convenience we summarize the singular 
moduli used. 



1 

= V2 - 1 
/C3=^V2(V3-1) 

k^ = {Vs-V2)i2-Vs) 

k^ - lV2{3-V7) 

o 

k^Q = {VW - 3){V2 - If 

k,, = -1=(2 - v^)(V5 - v^)(3 - v^). 
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and for K{kj^), correspondingly, 



4' 
93/4 

Kik^) = ^f3 (1/8) 
21/331/4 

21/1201/4 

= ,0 (V2-l)(V3+l)/3(l/24) 



48 

2^/773/4/? (1/7) /5 (2/7) 



14 /?(1/14) 



where the beta function, 



l3{x) = B{x,x)- 



T{2x) ' 
has been used. 

The energies can be calculated at the values C = Cjv = 1/VN- We give below 
the first three values, the simpler remaining ones being relegated to Appendix C. 
The scalar closed expressions in three dimensions are, from (18) and (45), 

r8(i/4) 



^3(^1) 
^3(^2) 
^3(^3) 



5120 7r6 

(3-2 v^) r^(l/8) 
3072 7r2 r4(5/8) 

r4(i/6)r^(i/3) 

12288 TT^ 



and in five dimensions, 
^5(^1) = 



61440 



r^d) ((7 - 5 V2) r^d) + 24 (-3 + 2 v/2) r^(i)) 

^^^^^ 7077888 7r6 re (i) 

_ (r^(i) (56 2lrr^ + lir^(i))) 

■^^^^3''" 24772608 7r9 
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The spinor energies are given by (44). In one dimension, 



^1(^3) = 



967rr2(|) 
967r3 



and in three 



4(^2) = 



2560 7r6 

r4(i) ((-21 + 16v/2) r4(i) + 32 (8-5v^) n'-r^{})) 



49152 TT^r^d) 



In five dimensions, simpler spinor examples are, 

f 32r^(|) 
^5(^1) — 

_f (336 23 7r8 + 70 2i7r^r6(i) + 5ri2(i)) 

^^^^^^ " 96768 7r6 re (|) " 

We now present some numbers computed from the above formulae which, in- 
cidentally, agree well with a direct evaluation using the elliptic numerical functions 
available in Mathematica. As a check, for increasing N, i.e. decreasing tempera- 
ture, the values approach the standard Casimir values, which are usually expressed 
in Bernoulli numbers. This could be thought of only as a necessary confirmation 
that the second terms in (8) and (20) must tend to zero as ^ becomes smaller. 
However, it does mean that however complicated the closed forms become for large 
A^, they must ultimately tend to a Bernoulli expression i.e. a simple rational num- 
ber! The zero temperature quantities are completely entangled in the finite form 
expression. 
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For scalar s, 



^3(^1) = 0.006065678717786288844 

Eg (^2) = 0.004305183178876425232 

E^{Q = 0.004185448373868254025 

^3(^6) = 0.004166873656512456887 

^3(^7) = 0.004166726978035004716 

E^{^^q) = 0.004166669015900687631 

E^(^^^) = 0.004166666693680499456 

^3(^00) = 0.00416 
= 1/240 

For spinors, 



^5(^1) = -0.0005054732264821907370 

^5(^2) = -0.0005125278116483441533 

^5(^3) = -0.0005125654321803054872 

^5(^6) = -0.0005125661374804481332 

^5(^7) = -0.0005125661375588626478 

^5(^10) = -0.0005125661375661265283 

^5(^15) = -0.0005125661375661375647 

E^{^^) = -0.000512566i37566 
= -31/6480. 



5(^1) 
5(^2) 
5(^3) 

5(^6) 

^3(^7) 

4 (^10) 

^3(^00) 



= 0.03639407230671773306 
= 0.03543620656643122772 
= 0.03541764326165557577 
= 0.03541666779673327755 
= 0.03541666684440446497 
= 0.03541666666803304403 
= 0.03541666666666835152 
= 0.035416 
= 17/48 



^1(^1) 
^'1(^2) 
^'1(^3) 
^1(^6) 

^i(eio) 

4(^15) 

^i(eoo) 



= 

= -0.06007307862716140024 
= -0.07470339906473023869 
= -0.08242382467179113065 
= -0.08242382467179113065 
= -0.08242382467179113065 
= -0.08332293842073511712 
= -0.083 
= -1/12. 



These numbers also show that the range, < ^ < 1, is essentially a Casimir 
range. Inversion can be used to obtain values in the complementary 'high temper- 
ature' region, 1 < ^ < 00. 



10. Conclusion and extensions. 



We have presented expressions for the internal energies of conformally coupled 
scalars and spinors on odd spheres in terms of elliptic function quantities. They can 
be used for convenient numerical evaluation, but our main aim, which has become 
somewhat overshadowed, was really to investigate the question of temperature in- 
version symmetry. Regarding this, it was only for the three-sphere that an overall 
symmetric quantity (the total energy) could be found. For the others only the 'par- 
tial energies' satisfied the symmetry. The same structure is noticed by Kutasov and 
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Larsen, [46]. Thus the statements in Cardy, [3], while correct, are, possibly, slightly 
misleading. See also Oshima [47]. 

The closed form expressions at temperatures corresponding to singular moduli 
are perhaps somewhat academic. They correspond to denominators in the thermal 
^-function, [1], of the form + N m? indicating number theory involvement, [48], 
and hence elliptic functions. Everything could be obtained from the Epstein form 
of the thermal ^-function. (See related calculations by Kennedy [49] on S^/F.) 
Indeed it was by the use of the known 'exact' forms of certain Epstein functions that 
Zucker derived values of K{kj^). Most other calculations also use Kronecker's limit 
formula, following Selberg and Chowla, [38]. It thus seems that our discussion is a 
little circular but it is systematic and has the advantage of uniformity of treatment. 
One must also not forget that there seem to be purely algebraic ways of obtaining 
K{kj^), [36], independent of the limit formula. Not only that but it would be 
otherwise difficult to treat K'/K of the more general Abel structure (50). 

As further justification of our elliptic approach, we draw attention to the fact 
that the internal energy on the arbitrary odd sphere, S'', is a polynomial in just two 
quantities, see (18) and (44). 

For scalars this is a simple consequence of standard analytical properties of the 
Weierstrass function, p, but seems not so obvious in the Epstein formulation. 

For spinors, the explicit construction of a cubic recursion seems to be new. It 
is most easily derived directly from the differential equation satisfied by ds. 

There are a large number of relations that we have not brought into the present 
discussion. For example, the double sum representations of thermodynamic quan- 
tities obtained in the earlier work [50] and [25] can be manipulated. One of the 
summations can be performed yielding a sum of inverse powers of hyperbolic func- 
tions. These constitute adequate numerical expressions and also can be turned into 
elliptic function form to make contact with the present approach, cf [51]. 

Regarding even spheres, no direct relation to elliptic functions could be de- 
tected. The sums relevant for standard thermodynamic scalars (bosons) are 




An obvious generating function, in the spirit of odd spheres, is 
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with periodicity 4iC, but no periodicity in the iK' direction. As it stands, this 
makes a connection with Jacobian eUiptic functions impossible and again indicates 
the necessity of altering the thermal periodicity. The same statements hold for 
spinors. In future work we intend to investigate this question. 

We will also investigate a general elliptic description of the complete thermo- 
dynamics of the system, not just the internal energy. 

Finally, regarding references, we have not been able to consult some classic 
texts, such as those by Enneper, Durege and Jordan on elliptic functions. 

Acknowledgment: KK was supported by the Max-Planck-Institute for Math- 
ematics in the Sciences, Leipzig. 
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Appendix A. Fermion expressions. 



As might be expected in a topic involving elliptic functions there are many 
routes to required results. This will be illustrated in this appendix by relating 
fermion and boson quantities. While not absolutely necessary, the analysis usefully 
invokes several basic pieces of elliptic technology. There are various ways of writing 
and re-arranging the fermion quantities. 

In terms of the Weierstrass approach via the twisted p, (27), one can seek to 
bring back the ordinary Weierstrass function, p, by re-arranging the summation 
lattice in (27). If this is drawn out, it is seen that the summations divide into a 
positive and a negative square lattice, displaced, one from the other, by, say, along 
the real axis and both rotated by 45° with respect to the original , 777,3 lattice. 
Therefore define new summation integers by M = rrii + 7713 and m = nii — 7773 • The 
positive lattice corresponds to M and 777 both being even, and the negative lattice 
to M and m both odd. Simple algebra then gives, 

p{u;u}-^,uj2) = p(w;a;i,a;2) - p(w + o^i; oii, u^a) + p{uj^;oj^,lo2) , (51) 
where _ 

iU-, — LUt — Wo 

(52) 

so that = (c!;2+5]^)/2 is a half-period of p{u; u^^^uj^). This function has a negative 
discriminant. The scaled rotation, (52), has determinant 2. 

The same transformation can be applied to the normal p-function, to give 

p{u] cc^i, UJ2) = piu; cDi, CD2) + Pi'^ + '^1; ^1,^2) - Pi^i,^!, ^2) ■ (53) 

Combined with (51) this yields, 

p{u] ^2) = 2p{u] Wi, 002) - p{u; c^i, ^2) , (54) 

which can be used to give an expansion of the fermion pu in terms of that for pu. 
Thus, similarly to (13), we write, 

^ 00 

p{u; ^1, ^2) = ^ + y](2t - 1) «2*-2 ^ (55) 
^ 7^2 

so that, 

^ 00 

p{u; oj^, u;2) = ^ + J](2t - 1) (2G, - G,) u^'-^ , (56) 

t=2 
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and therefore 



(57) 



allowing r]^ to be found from (23). This will be taken further in Appendix B. 

Relation (54) expresses the fermion function, p, in terms of the boson one, p, 
which is sometimes useful. 

Introducing yet another, more flexible notation, equivalent to the zs function 
is the theta function construction (c/ [19], vol. Ill p. 117), 

2Kzs{2Ku) = $(u,r) = log^i(w, r) , (58) 

du 



with 



Thomae denotes $ by Z^^, [27] p. 295. 

In terms of the fermion function, ds, can be written formally as the difference 
of two composite expressions, (c/ [13], §§84, 85), 

2Kds{2Ku) = ^$(^2, (1 + r)/2) - ^$((« + l)/2, (1 + t)/2) 

- d eM2Al + r)/2) 
du ^ ^2(^2,(1 + r)/2) 

_ d e,{u/2,T)e,{u/2,T) 

du ^ 6'2(w/2,t)6'o(w/2,t) ' 

The last equality is a consequence of the quadratic Gauss transformation expressed 
in theta functions, {e.g. [19] II, pp. 119, 120, 207). This can be confirmed by combin- 
ing the relevant qf-series by eye {e.g. [19] III,p.ll7, [20], p. 431, [18]) and this leads 
onto another way of relating fermion and boson expressions which is to start from 
the left-hand side of the relation (22), and work at the q-seiies level by applying, d 
la Jacobi, the trivial, and common, identity, 

X X 2a;2 ^^^^ 



1 + X 1 — X 1 — x-^ 
to the quantities r]f{^) of (20), which we can write as 



^ _ o2t-i ^2t(l/2) ^ i2n + ir-W^+' 
VA^)-^ 4t +^ l + g2n+i 

n=0 
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The sum over odds is written as, 

(2n + l)2^-lg2n+l ^ ^ ^2t-lgn {2n)^t-lq2n 

/ -I \ _|_ q2n+l Z-^ 1 ^ qn Z-^ 2 + Q^" ' 

n=0 n=l n=l 

and (60) used on both parts, to give for the right-hand side, 

2-^ \ — qU / J Y — Q^n / ^ ]^ _ (;2n / ^ ]^ _ (y4n 

n=l n=l n=l n=l 

Then, simple algebra, yields the relation, 

= ^t(20 - 2e,(a - 2^^-1(6,(0 - 26,(02)) 

= 0,(20-20,(0 (61) 

= 2,(0-2(1 + 2^^-2)6,(0, 

where 0,(0 = 6,(0 - 22*-ie,(02) and S(0 = e,(20 +22%(02). The combination 
(61) possesses the correct inversion behaviour, as a check. The two parts on the 
last line possess the inversion property separately. In terms of the right-hand 
side involves ^ q^l'^ and q. 

This expression can also be obtained by reorganising the summations in the 
Eisenstein series, H^^^ and using the basic equation, (23). 

In the lowest case of the circle, t = 1, it is easily confirmed that the zero mode 

element of the scalar energy typically cancels on the right-hand side to leave exactly 

— f 

the same combination for the total fermion energy, , 

E/(0 = :Ei(20 - (0 + . (62) 

Some relevant calculations on the relation between p and <p can be found in 
Tannery and Molk, [19] IV, pp. 37-40. In our notation they have, 

p{u\uj^,uj^) = p{u\uj^,uj^ — ^ ^ — (63) 

which also follows from (51), and leads to another form for p, 

p{u; u^, u^) = p{u; + ~^i^i^2 ~ . (64) 

We have not used this although it does yield the required expansions after, say, 
the introduction of two Jacobi functions, e.g. sn and dc . 
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It is possible to avoid the new modulus, k, by repeating the transformation to 
the periods, u}i,uj2- This effectively amounts to doubling the periods to 2u!2,2u!i, 
which are periods of p. The necessary analysis is given by Tannery and Molk [19] 
rV pp. 276-280. Their quantities y, Y and Z are, in our notation, 

Y = p{u;uj-^,uj2) , y = p{u;uj-^,uj2) , Z = 2a;i, 20)2) , 
so that the fermion p is. 

Now y can be expressed rationally in terms of Z which leads to, 

p^2Z + 2 ^^^~^'^^^^~^^^ -y, (65) 
Z - 63 

where the are the roots associated with Z. By homogeneity Cj = 6^/4 and clearly 
the modulus associated with Z remains k. Equation (65) can provide another way 
to expand the p. 

As a side remark, since p{u) has periods 2u)i, 2lj2, it can be expressed rationally 
in terms of Z. We find, using the expression for p{2u) in terms of p{u), 

~. 4^^ + 472^^ + 1073^ + 71/4 Je2-e,){e2-e^) 

p(u; 0)1,^2) = -^3 h2 

4Z^ - 72Z - 73 Z- 62 

where 72 and 73 are the invariants associated with Z. Homogeneity allows every- 
thing to be written in terms of p^ = p{u/2; lJi,lj2)- Thus 

^/ N 1 1 1 

and so 

-/ ,. , , , ^ 1 4p^ + 4g2pg + lOggPo + 92^4: (62 -ei) (62 -63) 

p[u,u;i,u;2) — 1 —, r . [Ob) 

4 4p3 _ g^p^ - 2(p„ - 62) 

This does not seem to have any computational advantages. 
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Appendix B. The power series expansion of elliptic functions. 



For our purposes in this paper, we need the expansions of zs and ds. Glaisher, 
for example, Hsts all the expansions but it might be useful to present an outline of 
some derivations and other pertinant information. 

The expansion of the elliptic functions sn u etc. as power series in u formed 
the subject of Weierstrass' first major work dating to 1840. The derivation is a 
technical issue that is not treated in much detail in most textbooks. Exceptions 
include Pricke, [52] and Briot and Bouquet [45]. 

There are many equivalent approaches depending on the representations of the 
elliptic functions employed. For example, one could use the theta form of zs , given 
later, or one could employ the relation between zs and the Weierstrass ^-function, 

. r]u , ri 
zs u — Cu — — , zs u — — (p u + — 
K K 

The calculations for sn, cn and dn are described in Fricke, [52] pp. 396+, and 
Konigsberger, [53] , following Weierstrass. There are two basic approaches. The first 
is a direct application of Taylor series. By successive differentiation, Jacobi showed 
that, for example, 

d^^sn u 



and 



= (4") + ai^^sn^ U + ... + ai,^)sn2- u) sn u 



= (4"^ + SaS'^^sn^ u + ... + {2v + l)a^:hn^^ u) cn«dn 



u 



where the a^^^ are easily found by recursion which shows that they are integer 
functions of the squared modulus, A;^. Then, using the values snO = 0, cnO = 1 
and dnO = 1, the Taylor series of sntt follows directly. A similar technique holds 
for cn and dn and for the other elliptic function combinations. Some coefficients 
are given in [19] IV pp.92,94. 

Obviously one can just use the series for the three basic functions, sn, cn and 
dn, to construct those for the others, say for ds = dn/sn, by further recursion or 
inversion. One could also proceed directly, as for sn, using ds'^ = (/c^ + ds^)(ds^ — 

A related technique, and the one referred to in the main body of this paper, is 
to use the derivative relations (we give these for the case we are interested in) 

d d ^ d 

— asu= — cs tt • ns tt , — cs tt = — ns u ■ asu , — ns tt = — ds u ■ csu ^ 
du du du 
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mutually connecting the trio of functions, ds, ns, cs. 

The conventional way of extracting the information in these relations consists 
of setting up three coupled recursion relations for the expansion coefficients of the 
three functions, e.g. Fricke [52], p. 400. Alternatively, further differentiation yields 
the differential equation for ds, 

ds"u = 2ds^u - (1 - 2/c^) dsu, 

and substitution of the expansion (32) yields, after a little algebra, the cubic recur- 
sion for J > 3, 



_ 4(2j-l)! ^ 
(^•-2)(2j + l)^ 



(2Z-l)!(2j-2Z-l)! 

(67) 



j-i-i 

E 



4:VlVnVj-l-n 



^ (2/-l)!(2n-l)!(2i-2Z-2n-l)!_ 

This is the fermion analogue of the boson recursion, (14), coming from p. 

In content, (67) is the same as the three coupled recursions, but is more compact 
and clearly shows that all coefficients are determined by just the ffist two, r]i and 
r/2. It is quickly checked that the results agree with those listed by Glaisher and 
leads to (43) and (49). It may not be the most efficient way of determining the 
highest polynomials, but it is formally instructive. 

For completeness we mention the second approach, adopted by Weierstrass, 
[54,55] , who writes, following Abel, the elliptic functions as ratios of power series by 
representing them as ratios, effectively of theta functions but more handily of the 
earlier Abel functions, Al^{u), which are just factors different from theta functions. 

One has 

AlJu) AUu) , AUu) 

AIq{u) Aloiu) ^^o(«) 

Weierstrass' next step is to Taylor expand the Al^{u) in a double power series in k'^ 
and u using the differential equations that they satisfy. The coefficients are then 
determined by recursion. The calculation is straightforward and well adapted to 
machine evaluation. The result is the series, 

Al„(u) = J2yr-^c^— , MM = X; (-i)".A„ 



(2m)!' ' ^(2m)! 

m=0 ^ ' m=0 ^ ' 
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with Aq = Bq = Dq = 1 and Cq = —1, = 0, from which the expansions for any 
eUiptic function can be found, after a further recursion in the form of a forward 
substitution. For example, writing, c/(32), 

one finds the relation between coefficients, 

£ ( 2z + 1 ) ^^-^ A = i'2r + 1)D,-A,, (69) 

and so ^ 

(2r + l)d, = (2r + 1)D, " " ^ (si + J) 

which is suitable for successive substitution. It is very easy to combine the recursions 
by machine, however Weierstrass [55] gives expressions for the A^, B^,C^, 
coefficients up to m = 10 and it is a simple matter to check that (70) reproduces 
the results from (67). 

This is an adequate scheme for practical purposes but somewhat utilitarian. An 
aesthetically more pleasing way of deriving the necessary series expansions is to use 
the equation for the twisted p-function, (56). In order to write this in terms of the 
modulus which corresponds to the original periods Ui,U2, the modulus, fc, asso- 
ciated with the rotated periods, cDj^, 0)2, has to be related to k, with a corresponding 
statement holding for the complete functions, K and K. The transformation from 
u)i, to cD^, is a combination of modular transformations with a Gauss quadratic 
one, c/ [19] II, p. 123, but a direct way of finding the relation is the following. 

Standard relations in elliptic function theory are 

k = J'-^, k' = J'-^ , — = v/i7^, 

V ei - 63 V ei - eg u-^ 

where the half-period values are, 

ei = P(^;^i,^2j ea = p( ^ \<^\,^2) eg = p( — ; c^i, c^a) . (71) 

and so the problem is reduced to finding the relation between the e^, and the e^,. 
This can be done from the connection (53) by substituting in various values of u 
and using the definitions (71). One has immediately that the last term in (53) is 
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just Putting u equal to uJi/2, lo^I'^ and a;3/2 in turn in (53), in the first two 
cases the <p functions on the right are evaluated at quarter periods. Using standard 
results, e.g. [20] p. 54, one quickly finds. 



Ci = 62 - 2iA/e2 - CgVei - 62 

63 = 62 + 2ii/e2 -63^61 -62 (72) 

^2 ~ ^1 "I" 63 ~ ^2 ~ ~2 62 , 

and it is then easy to show that the moduli and complete integrals are related by 

2kk' = -^, X = ikk'X. (73) 

It is to be remarked that the new modulus, k, is complex and that k' = k* so 
k"^ = 1/2 + ib with b real. 

It is straightforward to compute the fermion partial energies from (23) and 
(57). According to (57) one needs the combination 2g2 — §2 which, after using (73), 
gives, 

2fi^2-^2 = -y(8(A;A;')' + 7). (74) 

This is the same as in (49) , up to a factor. To deal with we write it in the form 
given by Ramanujan [23] eqn.(60), 

g, = Aa3(i + ^ikkr)Vl-{2kk'r. (75) 

This enables the transformations (73) to be simply applied in order to evaluate 
2g^ — ^3 which again reproduces the expression in (49). Note the useful symmetry 



1 - (2kk')'^ = a/1 - (2kk')'^ . 

2kk' ' 

From their structure, the quantities and rj^ will contain a polynomial in 2kk' 
for t even and a similar polynomial times ■>/! — {2kk'Y for t odd. This is actually 
very handy when evaluating at singular moduli because the combination, 2/c/c', is 
often the most convenient one to emerge from the modular equation, e.g. [33]. 

In terms of the expansion coefficients in (68), T inversion, (33), implies (com- 
pare with the scalar case, (48)), 

d^{k^) = {-IT d^{i - e) , 
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and so a more appropriate variable would be 1/2 — A;^ = k. For odd n, vanishes 
only at K = while if n is even, has a minimum at k = 0, maxima at k = ±Kq 
and passes through zero at k = ±K]^, where — > 1.5 as n — > oo. 

Incidentally, purely as a check with Glaisher, using the modular transforma- 
tions generated by T : k ^ k' and S : k' — > 1/k', the expansion coefficients for the 
functions cs and ns are given in terms of the d^. Those for cs are /c^"^ d^{l/k^) and 
those for ns, k'^"^ dj-k^ /k'^). This means that the d^ coefficients could also have 
been obtained by inverting the series for sn. 

Gudcrmann (Weierstrass' teacher) in 1839 produced expansions which used 
the transformation relations for the elliptic functions (see [27]). Hermite devised a 
related technique which is outlined in Briot and Bouquet [45], see also [19] II p. 209. 
One should also mention the work of Andre, [56]. More recently, a combinatorial 
formulation of the coefficients has been employed {e.g. Dumont [57]) which is more 
efficient for the highest coefficients. 
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Appendix C. Internal energies. 



In this appendix we display the remaining simpler singular case computations 
of the internal energies in three and five dimenions for scalars and spinors. This 
will give an indication of the complexity encountered, but the results are, at least, 
in finite terms. 

First for scalars, in three dimensions, 

_ (35- 20^2 + 16^3-14^6) (1/24) 

^^^^^^ " 138240 7r2 r4(13/24) 

i7r^(i/7) r^(2/7) r^(4/7) 

458752 ttS 

4 (7725 - 5460 ^2 + 3452 - 2442 ^/To) r4(9/40) r4(7/8) r4(41/40) 



^3(^15) 



7r2 r4(3/8) r4(21/40) r4(29/40) 
(245 - 104^5) r4(l/30)r8(4/15) 



8640000 23/5 ^2 r4(i/6) r4(l/3) r4(17/30) ' 
and in five, 

_ (3808 TT-^ r4(i) r^(f ) THI) + 171 ry^) r^(f ) r^d )) 

^^^^^ 1233125376 7ri2 

The corresponding spinor expressions are 



and 



f 2(l + 2v/2-2V3)r2(i) 

-^il^e) = 



e{{^i,) 



(r2(i)r2(f)r2(|)) 

128 TT^ 



(^32 -23 v^ + 4^ifp -90 V2) rHi)rHf,)T\i) 
i280 7r2r2(^)r4(|)r2(A) 
(4-v/5)r^(^)r4(^) 



36021 x/57rr2(i)r2(i)r2(|I) 



f ^ 224 r2(i) r^d ) r^d ) + 15 ry^) r^(f ) r^(|) 

^^^'^^ 114688718 
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